Persistence of Covalent Bonding in Liquid Silicon Probed by Inelastic X-ray Scattering 
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Metallic liquid silicon at 1787K is investigated using x-ray Compton scattering. An excellent 
agreement is found between the measurements and the corresponding Car-Parrinello molecular dy- 
namics simulations. Our results show persistence of covalent bonding in liquid silicon and provide 
support for the occurrence of theoretically predicted liquid-liquid phase transition in supercooled 
liquid states. The population of covalent bond pairs in liquid silicon is estimated to be 17% via 
a maximally-localized Wannier function analysis. Compton scattering is shown to be a sensitive 
probe of bonding effects in the liquid state. 
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Silicon (Si) presents a fascinating phase diagram as is 
the case in other systems that form tetrahedrally coor- 
dinated networks.pl Upon melting, Si transforms into a 
metal accompanied by a density increase of about 10%. 
The resistivity of liquid Si (Z-Si) at the melting tempera- 
ture T TO is 0.75 /ufim, which is comparable to that of sim- 
ple liquid metals such as Z-Al. However, the first neighbor 
atomic coordination number in Z-Si remains 5.5~6 [2j, 
which is approximately half that of simple liquid metals, 
hinting that covalent bonds survive even in the metal- 
lic state [3j. In fact, molecular dynamics simulations of 
molten Si at 1800K suggest that approximately 30% of 
the bonds are covalent and that these covalent bonds 
possess a highly dynamic nature, forming and breaking 
up rapidly on a time scale of 20 fs [4{. It is remarkable 
that two completely different types of bonds— metallic 
and covalent— can coexist in /-Si. In fact, the coexis- 
tence of two forms of liquid in a single component sub- 
stance has been predicted to undergo a phase transition 
as a function of temperature and / or pressure Q, and 
many theoretical studies support the existence of a liquid- 
liquid phase transition (LLPT) A recent study 
reports that Z-Si could undergo an LLPT below about 
1232K and above about -12kB, separating into a high- 
density metallic liquid (HDL) and a low-density semi- 
metallic liquid (LDL) Q. But, 1232K is far below the 
melting temperature of 1683K of Si, and as a result the 
supercooled state has remained inaccessible to current 
experimental techniques. Very recently, Beye et al. have 
performed time-resolved x-ray measurements on Si using 
a femtosecond pulse-laser Q to reveal liquid polymorphs 
of Si which could support an LLPT, but these experi- 



mental conditions are far from being ideal [6H8J so that 
the experimental confirmation of an LLPT in Si remains 
an open question. 

A key requirement for the possibility of an LLPT ob- 
viously is that the metallic and covalent bonds coexist in 
Z-Si. Although experimental investigations of the atomic 
configuration hint at the existence of covalent bonds in 
Z-Si, surprisingly, soft x-ray [13] and magnetic suscepti- 



bility measurements (11| of electronic properties so far do 
not support this viewpoint in that all four valence elec- 
trons in Z-Si appear to behave like free-electrons. Emis- 
sivity and thermal conductivity of Z-Si also exhibit a free- 
electron like temperature dependence It is clear 
thus that the existence of the covalent bonds in Z-Si is 
not well established. In this letter, Compton scattering 
experiments on Si are reported using a levitation tech- 
nique. Experimental results are interpreted in terms of 
full quantum mechanical simulations to show the persis- 
tence of the covalent bond in Z-Si. Moreover, by invok- 
ing a number of different first-principles approaches we 
show that Compton scattering spectroscopy is a sensitive 
probe of bonding effects even in the presence of disorder 
induced broadening of the spectrum in the liquid state. 

As background, we note that Compton scattering 
refers to the inelastic x-ray scattering process at large 
energy and momentum transfers. To date, the tech- 
nique has been applied primarily to investigate electronic 
structure and fermiology of solid-state systems, where 
the measured Compton profile is directly related to the 
momentum density distribution of the electronic ground 
state within the impulse approximation [l3l The 
present study shows how the Compton technique can also 
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provide a novel spectroscopic window on the liquid state. 
Since no charged particles entering or leaving the sam- 
ple are involved, Compton technique is a genuinely bulk 
probe, which is not complicated by surface effects present 
in photoemission or electron scattering experiments [l5| . 

Compton profiles of polycrystalline Si (300K) and 
molten Si (1787K) were measured by high-energy (116 
keV) inelastic x-ray scattering at the BL08W beam line 
of SPring-8 flol ] with a momentum resolution of 0.16 a.u. 
The data processing to deduce the Compton profile from 
the raw energy spectrum consists of the following pro- 
cedures: background subtraction and energy-dependent 
corrections for the Compton scattering cross-section, the 
absorption of incident and scattered x-rays in the sample, 
the efficiency of the analyzer and the detector, and cor- 
rection for double scattering events; see Supplementary 
Material for details of data handling procedures. Back- 
grounds of the measured profiles are suppressed in the 
levitation technique used since no material is present in 
the immediate vicinity of the sample at the x-ray scat- 
tering center. 

Molten Si is highly reactive with most crucibles. To 
hold the sample without contamination, a high temper- 
ature electrostatic levitator (HTESL) was used. The 
HTESL levitates a spheroid sample of 2mm diameter 
in a high-vacuum environment (approximately 10~ 5 Pa) 
usin g el ectrostatic forces via a feedback computer con- 



trol [17j, Ll8j. Two separate Be windows (1.5 mm thick) 



were installed in the HTESL chamber for the entrance 
and exit of the x-ray beam. The scattering angle was set 
at 165 degrees. The sample (Si, 99.9999% purity) was 
heated and melted using the focused radiation of three 
50 W semiconductor laser beams emitting at 808 nm. 
Temperature was controlled within 15K and measured 
by pyrometry with wavelengths of 0.90 and 0.96 nm. 

Car-Parrinello Molecular Dynamics (CPMD) simula- 
tions ;19] were performed with the Quantum-espresso 
package [2(| within the framework of the density- 
functional theory using the generalized gradient approxi- 
mation (2lj . We employed norm-conserving pseudopo- 
tentials with plane-wave expansion of the Kohn-Sham 
wavefunctions and charge density up to a kinetic energy 
cutoff of 24 Ry and 96 Ry, respectively. Spin-unpolarized 
simulations were done with 64 Si atoms in simple cubic 
supercells. A time step of 0.097 fs was used with the ficti- 
tious mass of electron set to 300 a.u., which is within the 



suggested range for accurate CPMD simulations \2'A [23j • 
The solid Si (s-Si) simulation was performed at 300 K at 
the experimental density of 2.33 g/cm 3 . The /-Si simula- 
tion was performed at 1787 K at the experimental den- 
sity of 2.56 g/cm 3 [H]. Nose-Hoover thermostat [l^, [26 1 
on ions was applied in both solid and liquid simulations. 
We also used an additional Nose-Hoover thermostat on 
electrons in the liquid simulation due to the metallic 
character of the system. The solid and liquid simula- 
tions lasted for around 9 ps and 13.6 ps, respectively, 
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FIG. 1: Experimental valence-electron Compton profile of 
solid (blue open circles) and liquid Si (red solid circles). The 
error bars are less than the symbol size. 



after thermalization. The s-Si Compton profile was cal- 
culated from the average of ten uncorrelated configura- 
tions, equally spaced by 1 ps. The Z-Si Compton profile 
was calculated from the average of ten uncorrelated con- 
figurations, equally spaced by 1.5 ps. Both profiles were 
spherically averaged and convoluted with a Gaussian to 
match the 0.16 a.u. experimental resolution. In order to 
estimate disorder effects in the liquid [27], we have per- 
formed the simulation on Z-Si using the first-principles 
Korringa-Kohn-Rostoker coherent potential approxima- 
tion (KKR-CPA) framework [28[. In these simulations, 
we adopted the lattice model for Z-Si given in Ref. [29j |. 
In particular, we used a body centered cubic structure 
with lattice constant 5.13 a.u. in which lattice sites are 
occupied by randomly distributed 50% Si atoms and 50% 
empty spheres. 

The valence-electron Compton profiles of s-Si (300K) 
and /-Si (1787K), obtained by subtracting the theoret- 
ical core-electron profile J c {Pz) from the experimental 
profiles, are presented in Fig. 1. The theoretical core- 
electron profiles are based on the free-atom Hartree-Fock 
simulations 30], where (ls) 2 (2s) 2 (2p) 6 are treated as core 
electrons. 

Figure 2 presents the differences between the Comp- 
ton profiles for the solid and liquid phases [AJ(p z ) — 
Jsoiid(Pz) - JiiquidiPz))- The results of CPMD simula- 
tion are obtained by taking the difference of the liquid 
Compton profile and the spherically averaged profile of 
the solid. The CPMD results agree very well with the ex- 
perimental results, establishing the efficacy of the bond- 
ing properties of /-Si obtained through CPMD simula- 
tion The difference between the KKR-CPA results 
for /-Si and the spherically averaged s-Si is also shown 
in Fig. 2. The effect of disorder in the KKR-CPA calcu- 
lation yields a momentum density smearing near p z = 1 
a.u. [28[ and produces a dip of A J(p z ) at this same loca- 
tion. However, in the experiment AJ(p z ) shows a peak 
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FIG. 2: Differences between the Compton profiles for the 
solid and liquid phases, [J so iid (j>z) - Jliquid (pz)]'- experiment 
(black), CPMD simulation (red), and KKR-CPA simulation 
(blue). 



near p z = 1 a.u. instead of a dip. We therefore conclude 
that the features induced by the Si site disorder are over- 
shadowed by competing effects due to the covalent bond 
breaking. 

In order to analyze the bonding character in l-Si 
obtained through CPMD simulation, the maximally- 
localized Wannier functions (MLWF) analysis method 
is used [3lj]. MLWF's are constructed through unitary 
transformation of Kohn-Sham wavefunctions such that 
the average spread (variance) of the Wannier functions is 
minimized. Same as the Kohn-Sham wavefunctions, each 
MLWF represents two electrons in spin-unpolarized cal- 
culations. In the s-Si simulation, the spreads of MLWF's 
are all smaller than 1.84 A 2 and all the corresponding 
Wannier centers (WC) are in the middle between two Si 
atoms, showing that all electron pairs in s-Si are cova- 
lent. On the other hand, MLWF's calculated from the 
Z-Si simulation exhibit diverse characteristics. In partic- 
ular, a majority of MLWF's have spreads larger than 1.84 
A 2 , up to over 7.8 A 2 . The spread is an indicator of the 
spatial extent of the MLWF. The inability of the MLWF 
procedure to localize some of the Wannier functions is 
expected as the system is metallic and the MLWF's with 
large spreads contribute to the metallic properties of l- 
Si. We refer to such electron pairs as diffuse pairs. The 
MLWF's with small spreads mostly correspond to the 
covalent electron pairs in l-Si. However, there is a spe- 
cial type of MLWF whose spread is comparable to that 
of the covalent electron pairs but the WC is not shared 
by any two Si atoms, so that the corresponding electron 
pair may be called a lone pair. The lone pair MLWF is 
computed to have a spread of 2.21 A 2 32 1. 

In order to categorize different types of electron pairs 
in l-Si, we sort the MLWF's into three groups: covalent 
bond pairs, lone pairs, and diffuse pairs. All electron 
pairs with Wannier spread larger than 2.21 A 2 are as- 



signed to be diffuse pairs. Among electron pairs with 
spreads smaller than 2.21 A 2 , an electron pair is a co- 
valent bond pair if the following two criteria are satis- 
fied: (1) The WC lies between two Si atoms separated 
by less than 3.1 A, which corresponds to the first min- 
imum of the Si-Si radial distribution function; (2) The 
WC lies within 1.36 A from the mid-point of the two Si 
atoms. The 1.36 A equals square root of 1.84 A 2 , which 
is the largest variance of covalent bond pair MLWF's in 
s-Si. With this scheme, we found that on the average, 
there are 17%, 83%, and less than 1% of covalent bond 
pairs, diffuse pairs and lone pairs, respectively (33j |. Fig- 
ure 3 shows a snapshot from the simulation carried out 
on l-Si. Note that the population of different electron 
pair types depends on the criteria adopted. The covalent 
bond pair population of 17% can be taken as the lower- 
bound. This is because the choice of the Wannier spread 
cutoff of 2.21 A 2 is too tight for a bond pair or a lone 
pair since MLWF's are expected to be more dispersed in 
the metallic liquid phase. 

The key point of our analysis is that the existence of 
covalent bonds in metallic l-Si is clearly confirmed, and 
that liquid Si is not homogeneous at the atomic scale. 
The coexistence of two different bonding natures in a 
single component liquid is a precondition that deeply su- 
percooled l-Si can undergo LLPT [f|. Our study thus 
support the possible occurrence of LLPT. 

The temperature of l-Si in our study is 1787K which 
is far above the predicted LLPT occurrence temperature 
of 1232K Q. To infer what kind of change in l-Si could 
occur and cause LLPT with decreasing temperature, it 
is useful to look at the temperature dependence of den- 
sity of l-Si. The density of l-Si has been measured be- 
tween 1350K and 1850K 24]; a quadratic behavior of the 
density as a function of temperature in the supercooled 
region has been observed. This is interpreted as the co- 
existence of covalent and metallic bonds [34| • Above the 
melting point, the temperature dependence of the density 
shows approximate linearity, indicating little heterogene- 
ity in bonding properties. In fact, our study shows that 
the majority of the electron pairs are in metallic bonds. 
In supercooled states, the ratio of covalent bonds could 
gradually increase with decreasing temperature since the 
temperature dependence of density deviates from linear- 
ity and becomes quadratic. The covalent bond pairs 
could assemble to form LDL domains. When temper- 
ature drops to 1232K where LLPT could take place, the 
domain size could exceed a critical size as in the case of 
nucleation in liquids, separating Z-Si into HDL and LDL. 

Finally, as stated in the introduction, there has been 
a longstanding controversy related to the interpretation 
of the nature of bonding in l-Si: atomic structure ex- 
periments have given indications of presence of covalent 
bonds 0, [35|, while other experiments suggest that all 
valence electrons behave like free-electrons 15- Our 
study resolves this controversy by establishing the per- 
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FIG. 3: Snapshot from the simulation on l-Si at 1787K ex- 
hibiting Si atoms (yellow), covalent bond pairs (green), lone 
pairs (red), and diffuse pairs (translucent blue). Note, bonds 
connecting Si atoms are only guides for the eye as they do not 
correspond to the chosen cutoff bond length of 3.1 A. 



sistence of a significant number of covalent bonds in Z-Si, 
even though a majority of valence electrons in /-Si exist 
as diffuse pairs. 

To summarize, we have carried out Compton scatter- 
ing measurements on s- and /-Si coupled with CPMD 
and KKR-CPA simulations. Comparison with the KKR- 
CPA results shows that Compton scattering remains a 
sensitive probe of bonding properties despite the smear- 
ing of the momentum density induced by disorder in the 
liquid state. The excellent agreement of the CPMD sim- 
ulation with the experimental results establishes the ex- 
istence of covalent bonds in metallic Z-Si. The contribu- 
tions of the covalent bond pairs, diffuse pairs, and lone 
pairs are estimated to be about 17 %, 83% and less than 
1%, respectively. Our study thus provides evidence for 
the coexistence of two distinct bonding species, metallic 
and covalent, which is a prerequisite for the occurrence 
of LLPT in deeply supercooled /-Si, and thus supports 
the possibility of such a phase transition. Future Comp- 
ton experiments with intense x-ray pulses will be needed 
to probe how the silicon bonding properties evolve from 
metallic to covalent character and how an LLPT may 
occur. 
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